A robust measurement protocol for Bohmian velocities 



F. L. Traversa,i'2,[*] Albareda,^ M. Di Ventra,^ and X. Oriols^^Q 

^ Departament d'Enginyeria Electrdnica, Universitat Autdnoma de Barcelona, 08193-Bellaterra (Barcelona), Spain 
'^Department of Physics, University of California, San Diego, La Jolla, California 92093, USA 

A protocol for measuring Bohmian velocities using two position measurements (each one modeled 
by a Positive Operator Valued Measure) separated by a finite time-interval is presented. The 
protocol is very accurate and robust so long as the first measurement uncertainty divided by the 
finite time-interval between measurements is much larger than the Bohmian velocity, and the system 
evolves under flat potential between measurements. The difference between the Bohmian velocity 
of the unperturbed state and the measured one is predicted to be much smaller than 1% in a large 
range of parameters. Counter-intuitively, the measured velocity is that at the final time and not a 
time-averaged value between measurements. 



The velocity of a classical object, requiring two po- 
sition measurements, is trivially implemented in many 
apparati which control our daily activity. On the con- 
trary, in the quantum world, the first position measure- 
ment implies a perturbation on the quantum system so 
that the knowledge of the velocity without perturbation 
is hardly accessible. One can minimize the back-action of 
the measurement on the system using the pre- and post- 
selected weak measurement, developed by Aharonov, Al- 
bert and Vaidman (AAV) [1 more than two decades ago 
and receiving increasing attention f?-;?! nowadays. As 
a relevant example, the spatial distribution of velocities 
of relativistic photons in a double slit scenario has been 
measured, and the associated quantum trajectories re- 
constructed [5]. Recently, Wiseman clarified the physi- 
cal meaning of such weak measurement of velocities in 
non-relativistic scenarios [3]. Using the weak AAV value 
[1], he showed that the velocity distribution constructed 
from an arbitrary pre-selected state and a post-selected 
position eigenstate with an infinitesimal temporal sepa- 
ration between measurements corresponds exactly to the 
Bohmian velocity of the unperturbed state [5] . 

Naively, one could think that such velocities distri- 
butions pushed orthodox quantum mechanics towards a 
contradiction because of the exact measurement of the 
position and momentum of an unperturbed state. Obvi- 
ously, such contradiction is not present because the men- 
tioned experimental and theoretical velocity distributions 
are computed from an ensemble of identical systems (not 
from an individual one). In other words, Wiseman's re- 
sults can be interpreted either as the orthodox hy- 
drodynamic velocity [3] or as a genuine measurement of 
the Bohmian velocity [10]. Both velocities are mathe- 
matically equivalent. Since these ensemble Bohmian ve- 
locities are free from fundamental controversies, one can 
envision their use (instead of the wave function) as a 
practical interpretative tool for studying general quan- 
tum phenomena [lll[12]- However, the practical condi- 
tions for measuring Bohmian velocities in a laboratory 
are different from the idealized scenario studied by Wise- 
man 3]. First, the weak measurement in a laboratory 
can be outside of the linear-response regime assumed in 
the AAV development [13j. Second, position measure- 



ments have a small but finite uncertainty, meaning that 
the post-selected state is not an exact position eigenstate. 
Identically, the time-separation between measurements 
must be finite. In this letter, using the Positive Operator 
Valued Measure (POVM) framework [121 to deal with 
two position measurements separated by a finite time, 
we extend the original Wiseman's conclusions about the 
measurable Bohmian velocities into practical laboratory 
conditions (free from previous idealized assumptions). 

The POVM appears as a natural modeling of a measur- 
ing process [11] when the laboratory is divided into the 
quantum system and the rest (including the measuring 
apparatus). Thus, the perturbation of the state due to 
the measurement of the first position be defined 

through the Gaussian measurement (Krauss) operators: 

W^^C^ J dxe'^""^^ \x) {x\ , (1) 

where is the experimental uncertainty. The measured 
position belongs to the set 371 of all possible mea- 
surement outputs of the apparatus. For simplicity, we 
assume 9Jt = M in a ID system, being the extension to 
the 3D spatial domain trivial. Then, the normalization 
coefficient Cw — (VvJ'O'tu)^^^^ is fixed by the condition 
/ dx^W^yV^ = /. Due to the unavoidable uncertainty 
on any position measurement, we consider an equivalent 
operator for the second position measurement of Xs'- 

dxe ^ \x) {x\ . (2) 

From a large set of measured positions, x^, at time t^ 
and Xs at ts — t^ + r, we construct the experimental 
velocity as: 

VeiXs,ts) = T~'^E[{Xs - X^)\Xs], (3) 

being E[{xs — Xuj)\xs] the ensemble average of the dis- 
tance Xs—Xw, conditioned to the fact that Xg is effectively 
measured. Since i?[a;s|a::s] = Xg, the theoretical compu- 
tation of the velocity does only requires evaluating 
E[a;u,|xs] using standard probability calculus, 

T-ir I 1 / dX-u^X'iuPiXii] n Xg) 

n^v^l^s] = (4) 
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with P{xw n Xs) the joint probability of the sequential 
measurements oixw and Xs, and P{xs) of Xg. Both prob- 
abilities can be computed from the Born rule, using ([T]) 
and ([2]), as: 

P{x^ n X,) = (*| Vi^t c/t^t5^t/^|y^ (5) 
P(a;s) = / n Xs). (6) 



being \^{t.^)) = j^*) the initial state. Strictly speak- 
ing, contrarily to the AAV expression [1], we are using 
a weak measurement without post-selection. The final 
state of the system (determined by the time-evolution of 
the initial state \^) and the measurement processes) has 
no relevant effect when computing ([s]) and ([6|. 

Let us now analyze P{xs) in detail by substituting Eq. 
^ and ^ into Eq. Then, we have 



n^s) = ci 



dxy^dx' dx" e ^ 



xmx'){x'\UlSlSsUr\x"){x"\^). (7) 

Integrating over Xw and using Eq. ([2]), we can rewrite Eq. 
as: 



2 / / dx'dx" m x') e 



dxe 



{x"\ *) X 
x'\ Ul \x) {x\ Ur \x")) . (8) 



For a particle of mass m that evolves under a fiat poten- 
tial during r, we can evaluate {x\ Ur \x') using [15j 



{x\ Ur \x') = {in(2hT/m)y^^^ e (2"-/'") 



(9) 



Substituting Eq. ^ into ([s]) and solving the integral be- 
tween parenthesis, we have 



P{Xs) 



dx'dx" e '^^4' e-(w)'(^'-^")' 



{^x'){x'\Ul\x;){x,\Ur\x"){x"\^). (10) 



One easily realizes that the probability in (10 1 can be 



computed as P{xs) = U^SlSsUl \^) when the fol- 
lowing limit is satisfied. 



T mcr. 



(11) 



Let us emphasize that this condition, includes Wiseman's 
result [3] as a particular case: cr^ — > 00, CTs — >■ and 
r — ^ 0. Our development will justify the effective mea- 
surement of the Bohmian velocity (up to a negligible er- 
ror) for a broad range of a^, Cs and r. 

Identical steps can be done for the evaluation of 
/ dxwXwP{xw n Xs) in Eq. Q. The only difference re- 
sides on the integration on Xw, which in this case gives 
(a;'-Fx")/2exp[-(a;'-x")V4cr2]. \]sing j dxx\x) {x\ = 



X, under the limit (11 1, we obtain J dx^x^Pixw H a;^) = 
Re((^'| UlSlSsUrX I*)) . Finally, we can rewrite Eq. ^ 
as: 



u 1 _ Romulsis^u^xm) 

r,[Xyj\Xs\ - (*|c/tsjs,(7,|*) • 



(12) 



Next, we define the following (averaged) position Xs = 

(^-l WrSlSsXUr / (*| WrSlSsUr |*), SO that USiug Eq. 

(12) and the commutator [Ur,x], we get: 



without any reference to W„,. To further develop Eq. 



(13), we evaluate the commutator [C/t,^] using the 
Maclaurin series for Ur- 



(14) 



where H = /2m -I- ^ is the system Hamiltonian with 
V a fiat potential at the spatial region where the wave 
function is different from zero during the time between 
measurements. No restriction on V for other regions 
and times. Given two operators A and B, it can be 
proven that [A",B] = ^^^^[A B]A"~3 _ Then, be- 

ing [H,x] = —ih/mp and [H,p\ — 0, the commutator 
[H" , x] gives: 



ihn 



pH 



n-l 



and substituting Eq. (15) into Eq. (14) we obtain: 



[Ur,x\ = pUr, 

ni 



(15) 



(16) 



without considering the limit r — 0. Using Eq. ( 16 ) and 
the definition ([2]), a straightforward calculation for the 
numerator of Eq. (flSl) gives: 



Re((*| UlSlSs[Ur,x\ I*)) = Tj{Xs.ts) 



tCI \ dxJ{x,ts)&M-{^s-xfla% (17) 



where J(a;, tg) is the standard quantum current probabil- 
ity density [11 . SimUarly, we define (^-j U\S\SsUr |*) = 
C2/dx|*(x,t,)|2exph(a;, -x)Va2] = |*(x„i,)|2 for 
the denominator. Finally, the velocity, defined as Eq. 
([l3| divided by r, gives: 



v(Xs,ts) 



(18) 



This expression is just the Gaussian-spatially-averaged 
current density J{xs, tg) inside a tube of diameter di- 
vided by the corresponding Gaussian-spatially-averaged 
probability \^{xs,ts)\'^ . 

Whether or not the Gaussian-spatially-averaged value 



(18) is identical to the Bohmian velocity depends on 
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the measuring apparatus resolution, i.e. as, and the de 
Broghe wavelength A associated to |^). Under the limit 



OS < A, 



(19) 



one can assume r) « ^'(xs, <s) for x G [xs — as,Xs + 
as], so that ^{xs,ts) ~ '^{xs,ts). Identically, J{xs,ts) « 
J{xs,ts) and Xs = Xs- Then, Eq. (181 directly recovers 
the Bohmian velocity v{xs, ts) ~ v with: 



v{Xs,ts 



J , ts ) 



(20) 




-0.5 n , , 0.5 
Position (|a,m) 



Let us mention that the consideration as ~ X and the 
momentum p = h/X implies h/(mas) ~ w in the limit 



(111. 



Let us now compute the velocity variance. Since Xs 
and r are fixed in Eq. (k5|, var{Vf,) = var{xw) /t"^ ■ Thus, 



var{x^ 



E[xl 



{E[xw\xs\)'^ where E[2 



de- 



fined in Eq. (|4| is obtained from Eq. ( 20 1 . The eval- 
uation of ^ dx.u,x'^P(xyj C\ Xs) follows identical steps as 
in the computation of P{xs), where again the only dif- 
ference resides in the integral in x^^ that now gives 
{al/2 + {x' + x"f/A) exp[-(a;' - x"Y/4,al. Using again 



J dxx \x) {x\ — X and / dx x" \x) {x\ 
suit, under the limit (11), is: 



nxiixs 



2 w 



1 Rc{{-i'\UlslSsU^x^\-i')) 

I 1 Ro({-9\xUlSlSsU^x\^}) 
2 (*|(7+S+S,(7^|*) 



which finally gives |16j 



.ar(.) = g + ^g^(.,) + 0(^) 



(21) 



(22) 



where Qb{xs) is the (loc al) Bohmian quantum potential 
^HH. Under the limits ([ll| and the term al/{2T^) 



in Eq. ( 22 ) will be orders of magnitude greater than the 
other two. For an experimentalist, this means that the 
presence of the quantum potential on the spatial fluctua- 



tions of Eq. ( 22 1 will be hardly accessible, and that var{v) 
provides basically the value a^ of the apparatus. Such 
variance can be used to evaluate the number N of mea- 



surements needed to obtain ( 20 ) with a given error e{N), 



using the well know result from the probability calculus 
that e[N) = ^Jvar{v)/\fN w a,J{T^/2N). 

In order to test how robust (i.e. how independent of 
ayj , as and r) is the possibility of measuring the Bohmian 
velocity in a laboratory, we compute the (local) error 
eyj{xs) = \ve{xs) — v{xs)\- The calculation is simply 
achieved expanding exp[— (x' — a;")^/4(T^] in Taylor se- 
ries by keeping the first three terms and repeating the 
above calculations for P{xs) and P(a:u, fla^s) (notice that 
the previous calculations are equivalent to an order zero 
Taylor expansion) obtaining [T^: 



Am?al, 



2(1 - Tdxv)dxP - rpdlv 



(23) 



FIG. 1. (Color online) Velocity distribution (n black solid 
line and «e red dashed line) and quantum current density 
(J black solid line and Je green dashed line) for an electron 
in a double slit experiment at ts — Hps, CTs = 0.2nm and 
= 150nm. Top inset, total error ea{xs) + Swi^s) for the 
velocity (black solid line) and current (green dashed line). 
Bottom inset, |vl>p at = Hps. 



where p = \i}{xs,ts)\'' 
apparatus error Ss{x. 
the requirement ( 19 ) 



We further define the measuring 
) = \v{xs) — v{xs)\ deriving from 
and a calculation similar to the 



the final re- previous one based on Taylor expansion gives [TB] : 



(Xs) 



"^d^p + {2d^p - pdx)dxv 



Ap + aldlp 



(24) 



It is worth noticing that, by construction, the total error 

£{xs) = \v{xs) — Vg^{xs)\ accomplishes e{xs) < £s{xs) + 

{Xs) • 

Finally, we observe that the same set of measured val- 
U.GS X and Xs can be used to define an experimental 
current density Je = t~^[P{xs)xs — J dxwXwP{xwf^Xs)]- 
We only have to change how the measured data is ana- 
lyzed. All the previous calculations can be then repeated 
for the current in a similar way. 

As a numerical test of our prediction, we consider an 
electron passing through a double slit. For simplicity, the 
time evolution of two ID initial Gaussian wave-packets 
with zero central momenta and central positions sep- 
arated a distance of 100 nm are explicitly simulated. 
This corresponds roughly to the evolution of the quan- 
tum state after crossing the double-slit at t = Gs. From 
Fig.[T] it is evident that the agreement between the ex- 
act Bohmian velocity v in ( |20| and the Vg [numerically 
evaluated from ([3]), (|4|, ([5]) and ^ without any limit or 
approximation] is excellent. In Fig. [2] we plot the nor- 
malized value of the error ew{xs) integrated over Xs as 
£u) — (/ dXs£w{xsY I J dxsv{xs)^y^^ ■ The main conclu- 
sion extracted from Fig.[2]is that a large set of parameters 
(large a^/T values) allows a very accurate measurement 
of the Bohmian velocity, justifying the robustness of our 
proposal. 

To conclude, we emphasize some relevant issues. First, 
we have shown theoretically and numerically that the 
Bohmian velocity of an unperturbed state under general 
laboratory conditions can be obtained from two POVM 
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(j„, (nm) 



FIG. 2. (Color online) Relative error e-w integrated over all 
positions Xs as a function of cr„ and r for as — 0.2nm for the 
numerical test represented in Fig.[l] Black line bounds the 
region for e{am) < 1% and red line is the analytical error for 
the value r = Ips. 
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FIG. 3. (Color online) Left inset, region of relative error 
Ew < 1% and right inset region of relative error Es < 1%. 
Solid lines are the boundaries for the velocity and dashed 
line are the boundaries for the quantum current. Dotted line 
bounds the a^/r region. 



measurements separated by a finite r. Unlike the re- 
sults derived from the AAV formulation [T], the limits 
(111 and (19) provide a simple quantitative explanation 



of the experimental conditions for an accurate and robust 
measurement of the Bohmian velocity. 



On the other hand, the error es{xs) in (24) has a term 
that diverges as (t'^/t, meaning that a t close to zero will 
produce an inaccurate measurement of the velocity for 
finite ■ This regime is reported in the right inset of Fig. 
[3j Roughly speaking, for t 0, the wave packet moves 
a distance vt. When vt < as the measured position Xg 
has no relation to the velocity. We emphasize again that 
Wiseman's result [3] does not suffer from this inaccuracy 
because he considers, both, CTs -> and r — > 0. 



A closer look at the expressions (23) and (24) shows 



that the error diverges when p has oscillations with min- 
ima tending to zero. This can be clearly seen in Fig.[l] 
where the highest peak of the velocity corresponds to a 
minimum of p very close to zero. This situation is re- 
versed when we evaluate the current J. In fact, in these 
critical points, J — > and even the corresponding errors 
become very small. In Fig.|3]it is evident the shift of the 
< 1% region due to this error reduction. 



Perhaps, the most surprising feature of our protocol is 
that a local (in time and position) Bohmian velocity can 
be measured with a large temporal separation between 
measurements, while one would expect a time-averaged 
value. This is highly counter-intuitive because we are 
in a scenario where the time-evolving interferences imply 
large acceleration of the Bohmian particle in order to 
rapidly avoid the nodes of the wave function. 

Finally, another relevant result is that the accuracy of 
the Bohmian velocity is obtained at the prize of increas- 
ing the dispersion on (as seen in Eq. ( 22 ) for large 
cr„). Therefore, the fact that we can obtain the Bohmian 
velocity is not because the system remains unperturbed 
after one position measurement, rather because of the 
ability of the ensemble average done in the Xw integrals 
on Eq. ([5]) and Eq. ([6| to compensate for the different 
perturbations. The fact that a very large perturbation 
of the state is fully compatible with a negligible error 
can be easily seen in our numerical data. The measured 
state is roughly equal to the product of the unperturbed 
wave function (whose support is i « 2000nm at time 
— lips in Fig. [I]) by a Gaussian function centered at 
the measured position with a dispersion equal to a^j (for 
example, cr^, w 150nm for r = Ips in Fig. [2]). Even for 
« L (i.e. a large perturbation), the velocity error is 
negligible in Fig.[2j 

The protocol presented here will be very welcomed by 
experimentalists because it clarifies the laboratory con- 
ditions necessary for an accurate measurement of the 
Bohmian velocity, while relaxing the experimental condi- 
tions needed for such measurements. As we have empha- 
sized in the introduction, the present work is fully com- 
patible with (and in fact, it is developed within) orthodox 
quantum mechanics. However, it opens relevant and un- 
explored possibilities of getting a deeper understanding 
of many quantum phenomena through the comparison 
between simulated and measured Bohmian trajectories, 
instead of using the wave function and its related param- 
eters. 
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